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We investigate the intermediate and late-time behaviour of the massive Dirac spinor field in the 
background of static spherically symmetric brane-world black hole solutions. The intermediate 
^ , asymptotic behaviour of the massive spinor field exhibits a dependence on the field's parameter 

D . mass as well as the multiple number of the wave mode. On the other hand, the late-time behaviour 

power law decay has a rate which is independent of those factors. 

CN ■ 

, PACS numbers: 04.50.-|-h 

^ \ I. INTRODUCTION 

Nowadays it is widely believed that extra dimensions play a significant role in the construction of a unified theory 
(~| \ of the four fundamental forces of nature. In such models it is often the case that our Universe can be treated as 
'i a submanifold to which the standard model is confined, embedded in a higher dimensional spacetime. If one takes 
^ ' the volume of the extra dimensions spacetime to be sufficiently large, one is able to lower the fundamental quantum 
, gravity scale to the electrovac scale of the order of a TeV. It is thus of interest to construct black hole solutions in 
' such brane-world models. The difficulties arising in such attempts stem from the fact that, in general, brane dynamics 
generates Weyl curvatures which in turn backreact on the brane dynamics. We can look at the problem in question 
CS| ' by projecting the Einstein equations onto the brane. This approach was introduced in Refs.[l|, It is also of interest 
CO I to think of a four-dimensional brane-world black hole solution as a slice that intersects a bulk black hole [1, 0, Q ■ In 
' Ref. ^ the possibility was raised of finding a regular Randall-Sundrum (RS) brane world on which a static spherically 
symmetric black hole, surrounded by realistic matter, is located. This was achieved by slicing a fixed five-dimensional 
OO , bulk black hole spacetime. On the other hand, studies of spherically symmetric brane-world solutions with induced 
■ gravity were extended to include nonlocal bulk effects Q . The scalar as well as the axial gravitational perturbations 
J> [ of what we shall call "brane-world black holes" were studied in Ref . Q . 

An important question for for black hole physics is the investigation of how various fields decay in the spacetime 
outside a collapsing body. The importance arises from the fact that, regardless of the details of the gravitational 
collapse and features of the collapsing body, the outcome of this process, i.e. the resultant black hole is characterized 
by just a few parameters such as mass, charge and angular momentum. The first researches in this direction were 
carried by Price in Ref.j9|] while the scalar perturbations on Rcissner-Nordstr0m (RN) background were considered 
in It was found that charged scalar hair decayed more slowly than neutral hair while the late-time 

tails in the gravitational collapse of a massive fields in the background of Schwarzschild solution were reported by 
Burko [l^ and in the the Reissner-Nordstr0m solution at intermediate late-time were considered in Ref. [16] . The 
very late-time tails of the massive scalar fields in the Schwarzschild and nearly extremal Reissner-Nordstr0m black 
holes were obtained in Refs. [l7j . (l8| . It was shown that the oscillatory tail of scalar field decays like t~^/^ at late 
time. The power-law tails in the evolution of a charged massless scalar field around the fixed background of a dilaton 
black hole were studied in Ref. [l^ , while the case of a massive scalar field was treated in (2^ . The analytical proof 
of the intermediate and late-time behaviour of the in the case of dilaton gravity with arbitrary coupling constant 
was provided in Ref.[2l|. On the other hand, the problem of the late-time behaviour of massive Dirac fields were 
studied respectively in the spacetime of Schwarzschild, Reissner-Nordstr0m and Kerr-Newman black hole [2^.[23l[2^. 
Ref. ^] was devoted to the analytical studies of the intermediate and late-time decay pattern of massive Dirac hair 
on a spherically symmetric dilaton black hole, in dilaton gravity theory with arbitrary coupling constant a. 

The growth of interests in unification scheme such as superstring/M-theory triggered in turn an interest in the 
decay of hair in the spacetimes of n-dimensional black holes. The no-hair and uniqueness property for static holes is 
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by now quite well established [26l|. The decay mechanism for massless scalar hair in the n-dimensional Schwarzschild 
spacetime was given in Ref. [27| . The decay pattern of scalar massive fields in the spacetime of rt-dimensional static 
charged black hole was discussed in Ref.[23|- It was shown that the intermediate asymptotic behaviour of the hair 
in question was of the form i-('+"/2-i/2) ^ Numerical experiment for n = 5 and n — 6 confirmed these results. In 
Ref. ^29,] the authors obtained fermion quasi- normal modes for massless Dirac fermion in the background of higher 
dimensional Schwarzschild black hole. 

As far as the brane-world black holes are concerned, Ref. ^39| was devoted to studies of the intermediate and late- 
time behaviour of the massive scalar field in the background of a static spherically symmetric brane-world black hole. 
Among other things, it was shown that the late-time power law decay rate is proportional to t~^^^. The massless 
fermion excitations on a tensional 3-brane embedded in six-dimensional spacetime were studied in [sij . 
The main aim of our paper will be to clarify what kind of mass-induced behaviour plays the dominant role in the 
asymptotic late-time tails as a result of decaying massive Dirac spinor hair in the background of brane-world black 
hole. 

The paper is organized as follows. In Sec. II we gave the analytic arguments concerning the decay of massive Dirac 
hair in the background of the considered black hole. Sec. Ill will be devoted to a summary and discussion. 



II. THE DECAY OF DIRAC HAIR IN THE BACKGROUND OF BLACK HOLE BRANE SOLUTION 



A. Spinor fields 



We shall begin our analysis by recalling the general properties of massive Dirac equation in an n-dimensional 
spherically symmetric background psj ] . Namely, we shall study the massive Dirac Eq. given by the relation 



7^V^V^-mj7/. = 0, (1) 

where is the covariant derivative = 9^ + \^lfjalb, M ^^'^ o, are tangent and spacetime indices. There are 
related by ej^, a basis of orthonormal one- forms. The quantity w"'' = are the associated connection one- forms 
satisfying de"' + wt" A e*" = 0. On the other hand, 7^ are Dirac matrices fulfilling relation {7°, 7''} — Irf-^ . 
If a metric takes the product form : 

g^.^dx^'dx'' = gab{x)dx''dx^ + g,nn{y)dy"'dy'' , (2) 

then Dirac operator _p satisfies a direct sum decomposition 

P=P^+P,. (3) 

If one defines a Weyl conformally rescaled metric by g^u = Vl^gi_iu one finds that 

pii) = V = (4) 

Because a spherically symmetric line element is necessarily conformally flat, a static metric, spherically symmetric 
metric of the form 



ds^ = -A^df + B'^dr'^ + C^dS^a (5) 

where A — A(r),B ~ B{r),C — C{r) are functions only of the radial variable r, and the transverse metric dYi'^^_2 
depends neither on t nor on r is conformal to an ultrastatic metric, one factor of which is conformally flat. This allows 
us to solve the Dirac equation by a succession of conformal transformations and direct sum decompositions. The 
assumption that ^' is a spinor eigenfunction on the [n — 2)-dimensional transverse manifold E, leads to the equation: 

p^^i = A*. (6) 

In case of {n — 2)-dimensional sphere the eigenvalues for spinor where found in Ref. [s^. They imply 




(7) 
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where 1 = 0,1,... 

Having in mind the properties given above, one may suppose that 

ptfj = rmP, (8) 

and take the form of the spinor ip to he : 

V'^TT^^X'^*- (9) 
If one carries out the exphcit calculations it turns out that : 

{l"dt+j'd,)x = A{m-^)x, (10) 



where we have denoted by 



dy = jdr, (11) 



the radial optical distance (i.e., the Regge- Wheeler radial coordinate). On the other hand, the gamma matrices 7°, 7^ 
satisfy the Clifford algebra in two spacetime dimensions. One should remark that having in mind a Yang-Mills gauge 
field Afj_, an identical result can be provided on the transverse manifold E. Namely, we have 

l>E,A.* = A*, (12) 

where a Dirac operator twisted by the the connection A^. 

Finally, if we take into account that ip has the form as cx e~*"* one obtains a second order equation for x, that is 



d'^X ,2 .2/ ^ 



+ uj\ = A\m--Yx. (13) 



B. Dadhich-Maartens-Papadopoulous-Rezania (DMPR) brane- world black hole solution 

We treat first the case of the static spherically symmetric black hole localized on a three-brane in five-dimensional 
gravity in Randall-Sundrum model [33|. Having in mind the effective field equations on the brane one gets the 
following brane- world black hole metric [ij: 

ds^ = -U-^^ -T^V + T (14) 



M^r M2r^^ fi ^ 



where g is a dimensionless tidal parameter arising from the projection onto the brane of the gravitational field in the 
bulk, Alp is a fundamental five-dimcnsional Planck mass while Mp is the effective Planck mass in the brane world. 
Typically, one has Mp Mp. In what follows we shall concentrate on the negative tidal charge which is claimed [l| 
to be the more natural case. Thus, the roots of goo = are respectively r_|- and r_. Namely, they imply 



Mp2 V V Af2 M2 



Expressing the negative charge as Q, for simplicity, we can rewrite the roots as follows: 



r± = M(^l±yi + ^j. (16) 

Our main aim will be to analyze the time evolution of a massive Dirac spinor field in the background of brane- 
world black hole by means of the spectral decomposition method. In Refs.[la],[3l] it was argued that the asymptotic 
massive tail is due to the existence of a branch cut placed along the interval —m < oj < m. Thus, an oscillatory 
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inverse power-law behaviour of the massive spinor field arises from the integral of Green function G{y,y' around 
the branch cut. Consider, next, the time evolution of the massive Dirac spinor field provide by the relation 



X{y-,t) = / dy' 



G{y,y';t)xt{y',0) + Gt{y,y';t)x{y',0) 



for t > 0, where the Green's function G{y,y';t) implies 
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9*2 dy^ 



V 



G[y,y'-t) = 5[t)5[y-y'). 



(17) 



(18) 



By means of the Fourier transform ,34] G{y,y';uj) = dt G{y,y' ;t)e^'^*' , Eq. lfTS)) can be reduced to an ordinary 
differential equation. The Fourier transform is well defined for Im lu > 0, while the corresponding inverse transform 
yields 



G{y,y';t) = 



2tt 



oo+ze 



duj G{y,y';uj)e~ 



(19) 



— oo+ie 



for some positive number e. By virtue of the above the Fourier component of the Green's function G{y,y']U!) can be 
rewritten in terms of two linearly independent solutions for homogeneous equation. Namely, it reduces to 



d^ 
dy 



2 +uj'-V]x. = 0, 



(20) 



where V — A ym - ^ 

The boundary conditions for Xi are described by purely ingoing waves crossing the outer horizon of the static 
black hole xi — e"*'^^ as y ^ —00, while X2 should be damped exponentially at Thus, X2 — e~^™ ~^ ^ at 
y-^cG. 

Suppose now that the observer and the initial data are situated far away from the considered brane black hole. Let 
us rewrite Eq. (|20p in the more convenient form using the change of variables 



1 - 



1/2 



1/2 ' 



(21) 



where i = 1,2. Then, we expand Ea. (j20p as a power series of r±/r neglecting terms of order C'((a;/r)2) and higher. 
We obtain the the following: 



dr2 



2 2 2aj2(r_|_ + r_) — m2(r+ + r_) + 2Am(l + r+) 



(22) 



A2 — 2Amr_ + m^r^r^ 



The solution of equation ([221) may be obtained in terms of Whittakcr functions. Two basic solutions are needed to 
construct the Green function, with the condition that | w |> m, i.e., xi — Ms,p^(2Lur) and X2 = Ws.fi{2Lur). The 
parameters of them imply 



fl ~ ■\/l/4 + A2 — 2Amr_ + m?r+r- , 

^ _ w2(r+ +r_) + Am(l + r+) - i^(r+ +r_) 



(23) 



- 2 2 2 



Having all this in mind, we reach to the following form of the spectral Green function: 



Gc{x,y;t) = 



1 

2^ 



dw 



Xi(a;,a)e") X2(2/, t^e""') Xi{x,u}) X2{y,Co) 



W(uj) 



(24) 
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where W{u}) stands for Wronskian. 

Let us analyze first, the intermediate asymptotic behaviour of the massive spinor field with the range of parameters 
M ^ r ^ i <C M/{mMY. The intermediate asymptotic contribution to the Green function integral gives the 
frequency equal io Co — 0{^J m/t), which in turns implies that 5 Using the fact that 5 results from the 1/r term 
in the massive spinor field equation of motion, it illustrates the effect of backscattering off the spacetime curvature. 
In the case under consideration the backscattering is negligible. Thus, we find the following: 



Mr(2M) r(i - A) 



1 gCSM+i)""* 



(rr') 



(25) 



where one applied the fact that wr <C 1. We also have in mind that /(w) can be approximated using the fact that 
Af (a, b, z) = \ as z tends to zero. Consequently, the resulting spectral Green function reduces to the form as 



G (r r'-t) = 23^-^r(-2A) r(i + A)r(A + 1) 



{rr') 



(26) 



Taking into account the limit when t 3> 1/m we conclude that the spectral Green function yields 



Gcir,r';t) 



23A-ir(-2/2) r(i + A)r(A + i) 
AV^ fiT{2fi) r(i - m) 



1 + e(2'i+i)'^M (rr')^+'^ m'^ r'^-'' cos(mi - |(/1 + 1)). (27) 



We remark that Ea. (j27|) exhibits an oscillatory inverse power-law behaviour. In our case the intermediate times of 
the power-law tail depends only on fl which in turn is a function of the multiple number of the wave modes. 

The other pattern of decay of massive spinor Dirac hair is expected when 5 3> 1, for the late-time behaviour. 
Namely, when the backscattering off the curvature is taken into account. Under the assumption that (5 3> 1, /('^) 
may be rewritten in the form as 



r(i + 2^) r(i-2A) , ,,1 
^ 2A ^'' ^ 



J2f,iV8dQr) J_2m(V8^) " h-^iVSd^) I^2f.{V86Q?) 



(28) 



(r(i + 2/2))2 r(-2A) r(i + M-^) 
2Ar(2A) ra-fi-s) 



(rr')2 K 



-2/1 



+ e(2'i+i)/2^(V8J^) /2A(V8fcr') 



where we used the limit Ms.fi{2ujr) w r(l + 2p){2Cjr)i (5~^ J/i(V85(I;r). One should notice that the first part of 
Eq. ([28|) ■ the late time tail, is proportional to t^^. It occurs that we shall concentrate on the second term of the 
right-hand side of Eq. ([28|) . For the case under consideration it can be brought to the form: 



(29) 



where the phase is defined by the relation 



while M is given by: 



1 + (-l)2Ae-2'rj<5 

2 ^ (^_^^2fig2-!TiS ' 



M = (r(i + 2A))^ r(-2A) 

2fi r(2^) ^ ' 



J2^(V8fc^) J2A(V8(5(ir')+/2A(\/8fc^) I2^.{V8SC^r') 



(30) 



(31) 



The saddle-point integration allows one to find accurately the asymptotic behaviour. This method is applicable in our 
case because of the fact that at very late time both terms e™* and e^'^'^ are rapidly oscillating, which in turns means 
that the spinor waves are mixed states consisting of the states with multipole phases backscattered by spacetime 
curvature. Most of them cancel with each others which have the inverse phase. The saddle-point is found to exist at 
the following value: 



ao 



TT ([j2 (r+ +r-) + Am(l + r-) - ^{r+ + r_) 
V2m 



(32) 
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Then, the resultant form of the spectral Green function yields 

{mt)^« sin(mt) x(r, m) x(r', m). (33) 

The above form of the spectral Green function concludes our investigations of the late-time behaviour of massive 
Dirac hair in the background of DMPR brane-world black hole. The form of it envisages the fact that the late-time 
behaviour of the fields in question is independent on the field parameter mass as well as the number of the wave mode. 
The late-time pattern of decay is proportional to —5/6. 



2V2 

71" 



.2/3 



2m^(r+ +r_) + 2Am(l + r+) -m^(r+ +r_) 



C. Casadio-Fabbri-Mazzacurati (CFM) brane black hole solution 



Our next task will be to consider a general class of spherically symmetric static solution to fivc-dimensional equations 
of motion by considering the general form of the line clement provide by the metric 



ds^ = -A{r)dt^ 



1 



B{r) 



^dn^ 



(34) 



Casodio et al. [3] obtained two types of analytic solutions by fixing either A{r) or B{r). The solution will be given 
in terms of the ADM mass M and the parametrized post-Newtonian (PPN) parameter j3 which affects the perihelion 
shift and the Nordtvedt effect [3^ . The momentum constraints are identically satisfied by the metric coefficients and 

the Hamiltonian constraints can be written out Setting A{r) = ( 1 — j the resulting metric yields 



ds^ 



1 ]dt^ 

r 



1 _ 3M 



^ _ 2M 



-dr^ +r^dn^, 



1 _ iM 1 

^ 2r I 



(35) 



where 7 = 4/3—1. A convenient form of the equation of motion for a massive Dirac field can be obtained by the 
transformation : 



Xi 



ML] 

2r I 



I — Ml 



1/2 



1 _ iM 

^ 2r 



1/4?' 



(36) 



where i — 1,2. As in the preceding section, let us expand Eq. (fT8|) as a power series of Al/r neglecting terms of order 
0{{uj/r)^) and higher. It then follows directly that one has 



dr'^ ^ 



2 2 

u! — m + 



tj2 a - 6 + 2Xm + A M m (3 - 7) - H-P 7 



(37) 



where a = (54-7) and b 



(7-3). 



Thus, the two basic solutions which are needed to construct the Green function, with the condition that | w |> m 
are given by Xi — Ms^fi{2Lur) and X2 ~ W^(5,/i(2'^f'), with the following parameters: 



jl= W1/4 + A2 + A Af m (3-7)--M2 j 



6^ 



ui^ a~ w? h + 2\m 
2^ 



~2 2 2 

Lo = m — uj . 



(38) 



The preceding section arguments can be repeated. The conclusion is that the spectral Green function of the interme- 
diate late-time behaviour of massive Dirac spinor fields with the new parameters of the Whittaker functions given by 
the relation ([38|) . Consequently, the next step will be to calculate the late-time behaviour of the considered field. It 
can be verified that the stationarity of the integral will be achieved for the parameter 



TT (cl|2 a — vn? h + 2Am) 
2\/2m 



(39) 
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By virtue of saddle point method, on evaluating the adequate expressions, we find that the spectral Green function 
provides the following: 



4Mm2 + 2Ar 



(mt) 6 sin(mt) ■0(r, m) ^{r\ra)^ 



(40) 



One can observe that the dominant role in the late-time behaviour is played by the term proportional to —5/6. 
On the other hand, let us consider that B{r) = ^1 — ^^^^ for the other model of brane-world black hole. It implies 
the following line element: 



7 



7-1 + Wl 



27MY 



dt' + 7 7 



(41) 



Next, let us change coordinates as follows: 

Xi = 



72^ 



7_l + ^/l_22M fi_22M 



(42) 



where i = 1,2, Then, expand Eq.pB]) as a power series of M/r neglecting terms of order 0{{uj/r)'^) and higher. It 
yields 



■ . . n o 47M(tj2 (1 + p2^ _ _^ 2Xm ~ 8mMXj - 4M2m272 
tj 7 p - m H 5 



€ = 0, 



(43) 



where p'^ = — 1)^ + 3. 

Ea. (|15)) can be brought to the form of Whittaker's equation. Two basic solutions are needed to construct the 
Green function. The additional requirement that | lu |> m, implies that they are of the form — Ms_fj_(2ujr) and 
X2 = Ws^fi{2ojr). The parameters of the Whittaker functions are given by 

/i= V1/4 + A2 -8mA/A7-4Af2TO2^2^ — [ — ^ ^ ^ ~ lu^ Y p^ . (44) 



2^ 



On the other hand, from the considerations presented in the preceding case, the stationarity of 2TrS ~ tot can be 
obtained for the parameter equal to 



ao 



TT {AjM{uj'^{l + p2) - m2) + 2Am) 



2\/2m 

Summing it all up, one obtains the asymptotic late-time spectral Green function in the form 
2V2 



(45) 



Gc{r,r';t) 



V372[(7_i)2 + 3] 



to2/3 (^)I 



47Mm2((7 - 1)2 + 3) + 2Ato 



(mt) 6 sin(TOi) x(r, m) x(^', 'Ti), (46) 



As in the previous cases the dominant role in the asymptotic late-time decay of massive Dirac hair in the spacetime 
of CFM brane black hole plays the oscillatory tail with the decay rate proportional to t^^^^. 



III. CONCLUSIONS 



In our paper we have considered the problem of the asymptotic tail behaviour of massive Dirac hair in the spacetime 
of various brane-world black hole solutions. Our main aim was to reveal what type of mass-induced behaviours play 
the main role in the asymptotic intermediate and late-time decay pattern of black hole hair in question. In our 
considerations we took into account two brane-world black hole solutions given in Refs.[ll, It was shown that in 
the case of intermediate asymptotic behaviour one gets the oscillatory power-law dependence which varied with the 
multiple number of the wave mode I as well as with the mass of the Dirac fields. As in the case of ordinary Einstein 
static spherically symmetric black hole spacetimes this pattern of decay is not the final one. At very late-times the 
resonance backscattering off the spacetime curvature emerges, which in turn is independent on angular momentum 
parameter and the field parameter m. The late-time asymptotic pattern of the decay is of the form t~^/^. It should be 
interesting to find a general proof of this pattern of decay for massive Dirac fields in the spacetime of static spherically 
symmetric black object. The investigations in this direction is in progress and will be published elsewhere. 
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